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Description 

This invention relates generally to digital x-ray imaging and, more particularly, relates to methods and apparatuses 
for reducing scatter in two-dimensional x-ray imaging and two-dimensional dual-energy x-ray imaging. 

Recent advances in semiconductor industry have resulted in the ability of fabricating large-format two-dimensional 
integrated detector arrays for x-ray detection. These arrays have on the order of one million detector cells and provide 
instant acquisition of two-dimensional x-ray images with exceedingly high quality. 

Scatter, which results from those x-rays that strike objects and deflect in random directions, has been a difficult 
and on-going problem in x-ray imaging using two-dimensional detectors. For example, in projection chest radiography, 
scatter typically accounts for between approximately 30% and 50% of the total amount of x-rays detected. 

Two-dimensional detectors are exposed to wide-angle random scatter. Scatter tends to degrade image quality, 
and to block the way for quantitative imaging. For example, with scatter present, three-dimensional image reconstruc- 
tion cannot be implemented. 

Embodiments of the present invention use the dual-energy x-ray imaging method as a foundation for the method 
of eliminating scatter- However, the dual-energy x-ray imaging methods of prior art are not suitable for this purpose. 

All the dual-energy x-ray data decomposition methods of the prior art have a common approach: linearization. 

It is a consensus that none of current linearization methods is suitable for use with million-pixel two-dimensional 
detector arrays. One object of the invention is to provide a substantially improved dual-energy x-ray imaging method 
suitable for two-dimensional detector arrays. 

Currently, there are two basic methods for reducing scatter in two-dimensional x-ray imaging. The first method 
uses an anti-scatter grid or improved grid devices to slightly relieve scatter effects on images. By increasing the air 
gap between the subject and the detector, the scatter can also be reduced, but at the same time, the image is blurred 
due to the geometric distance the x-rays have to travel. 

The second method to reduce scatter is to calculate theoretical estimates. Theoretical methods, including Monte 
Carlo simulation methods and analytical deconvolution methods, can only give crude predictions. 

According to a first aspect of the present invention, there is provided a two-dimensional x-ray imaging system for 
taking images of a subject, said system comprising: 

(a) an x-ray source adapted to emit x-rays for passage through said subject, said x-rays including primary x-rays 
having their direction of travel unaltered by interaction with said subject and said x-rays including scatter x-rays 
having their direction of travel altered by interaction with said subject; 

(b) x-ray detection means for detecting said x-rays in two dimensions; 

(c) said subject being located between said x-ray source and said x-ray detection means; 

(d) said detection means including high-resolution detection means for detecting said primary x-rays and said 
scatter x-rays simultaneously; 

(e) beam selection means for blocking passage of substantially all of said scatter x-rays and permitting passage 
of a portion of said primary x-rays; 

(f) said detection means including low-resolution detection means for detecting only said x-rays passed by said 
beam selection means. 

According to a second aspect of the present invention, there is provided a method for taking a two-dimensional x- 
ray image of a subject using a two-dimensional x-ray imaging system, said subject being composed substantially of 
two materials that interact differently with x-rays, said imaging system including an x-ray source, beam selection means, 
and a two-dimensional x-ray detector assembly having a plurality of detection locations, said subject being between 
said x-ray source and said detector assembly, said x-ray source being adapted to emit x-rays for passage through said 
subject, said x-rays including primary x-rays having their direction of travel unaltered by interaction with said subject, 
said x-rays including scatter x-rays having their direction of travel altered by interaction with said subject, said beam 
selection means having a blocking position wherein passage of substantially all of said scatter x-rays to said detector 
assembly is blocked and passage of said primary x-rays to said detector assembly is permitted, said beam selection 
means having a transmissive position wherein passage of substantially all of said primary x-rays and said scatter x- 
rays to said detector assembly is permitted, and selected detection locations being those of said detection locations 
receiving only said primary x-rays when said beam selection means is in said blocking position, said method comprising: 

(a) illuminating said subject with said x-rays; 

(b) moving said beam selection means to said transmissive position; 

(c) acquiring a raw high-resolution image l m from said detection locations and processing said image l^ to normalize 
it and to subtract dark signals, yielding a high -resolution image composed of said primary x-rays and said 
scatter x-rays; 
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(d) producing, from said image Dft,, a low-resolution image D fl representing said selected detection locations; 

(e) moving said beam selection means to said blocking position; 

(f) acquiring a raw low-resolution image l rt from said selected detection locations and processing said image | rt to 
normalize it and to subtract dark signals, yielding a low-resolution image D ri composed of substantially only said 
primary x-rays; 

(g) calculating a low-resolution primary x-ray image D fPI by multiplying said image D ri by a predetermined constant 
that is independent of said subject; 

(h) calculating a low-resolution scatter x-ray image D e , by subtracting said image D fP , from said image D fl ; 

(i) calculating a high-resolution scatter x-ray image D,s h by extending said low-resolution scatter x-ray image D^, 
to all of said detection locations through interpolation; and 

(j) calculating a high-resolution primary x-ray image Dfp h by subtracting said high-resolution scatter x-ray image 
Dfs h from said high-resolution image D^; 

(k) whereby said image Djp h is a high-resolution, two-dimensional primary x-ray image of said subject after said 
scatter x-rays have been substantially eliminated, said image Dfp h having a resolution substantially equal to the 
highest resolution available from said detector assembly. 

According to a third aspect of the present invention, there is provided a method for taking a two-dimensional x-ray 
image of a subject using a two-dimensional x-ray imaging system, said subject being composed substantially of two 
materials that interact differently with x-rays, said imaging system including, in physical sequence from front to back, 
an x-ray source, a front two-dimensional x-ray detector assembly having a plurality of front detection locations, beam 
selection means, and a rear two-dimensional x-ray detector assembly having a plurality of selected rear detection 
locations, said subject being between said x-ray source and said front detector assembly, said x-ray source being 
adapted to emit x-rays for passage through said subject, said x-rays having an energy spectrum that can be approxi- 
mated as having a single energy E 0 , said x-rays including primary x-rays having their direction of travel unaltered by 
interaction with said subject and said x-rays including scatter x-rays having their direction of travel altered by interaction 
with said subject, said beam selection means blocking the passage of substantially all of said scatter x-rays to said 
selected rear detection locations and permitting the passage of said primary x-rays to said selected rear detection 
locations, said selected rear detection locations receiving only primary x-rays, and selected front detection locations 
being those of said front detection locations intersected by x-ray projection lines extending from said x-ray source to 
said selected rear detection locations, said method comprising: 

(a) illuminating said subject with said x-rays; 

(b) acquiring a raw high-resolution image l m from said front detection locations and processing said image l m to 
normalize it and to subtract dark signals, yielding a high-resolution image composed of primary x-rays and 
scatter x-rays; 

(c) producing, from said image D^, a low-resolution image D fl representing said selected front detection locations; 

(d) acquiring a raw low-resolution image l ri from said selected rear detection locations and processing said image 
l ri to normalize it and to subtract dark signals, yielding a low-resolution image D rt composed of substantially only 
said primary x-rays; 

(e) calculating a low-resolution primary x-ray image D fP , at said front detector assembly by multiplying said image 
D rt by a predetermined constant that is independent of said subject; 

(f) calculating a low-resolution scatter x-ray image D fe , at said front detector assembly by subtracting said image 
D fP | from said image D fl ; 

(g) calculating a high-resolution scatter x-ray image D,g h by extending said low-resolution scatter x-ray image D e , 
to all of said front detection locations through interpolation; and 

(h) calculating a high -resolution primary x-ray image Dfp h by subtracting said high-resolution scatter x-ray image 

from said high-resolution image D^; 

(i) whereby said image Dfp h is a high -resolution, two-dimensional primary x-ray image of said subject after said 
scatter x-rays have been substantially eliminated, said image D fPh having a resolution substantially equal to the 
highest resolution available from said front detector assembly. 

According to a fourth aspect of the present invention, there is provided a method for taking a two-dimensional x- 
ray image of a subject using a two-dimensional x-ray imaging system, said subject being composed substantially of 
two materials that interact differently with x-rays, said imaging system including, in physical sequence from front to 
back, an x-ray source, a front two-dimensional x-ray detector assembly having a plurality of front detection locations, 
beam selection means, and a rear two-dimensional x-ray detector assembly having a plu rality of selected rear detection 
locations, said subject being between said x-ray source and said front detector assembly, said x-ray source being 
adapted to emit x-rays at two different average energy levels for passage through said subject, said x-rays including 
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primary x-rays having their direction of travel unaltered by interaction with said subject and said x-rays including scatter 
x-rays having their direction of travel altered by interaction with said subject, said beam selection means blocking the 
passage of substantially all of said scatter x-rays to said selected rear detection locations and permitting the passage 
of said primary x-rays to said selected rear detection locations, said selected rear detection locations receiving only 
5 primary x-rays, and selected front detection locations being those of said front detection locations intersected by x-ray 
projection lines extending from said x-ray source to said selected rear detection locations, said method comprising: 

(a) illuminating said subject with x-rays of a first average energy level H; 

(b) acquiring a raw high-resolution image \ Wh from said front detection locations and processing said image l^ 
10 to normalize it and to subtract dark signals, yielding an image which is composed of primary x-rays and scatter 

x-rays; 

(c) producing, from said image D^, a low-resolution image D^, representing said selected front detection loca- 
tions; 

(d) acquiring a raw low-resolution image l^, from said selected rear detection locations and processing said image 
is IfHj to normalize it and to subtract dark signals, yielding an image D^, composed of substantially only said primary 

x-rays; 

(e) illuminating said subject with x-rays of a second average energy level L; 

(f) acquiring a raw high-resolution image ^ from said front detection locations and processing said image l^ to 
normalize it and to subtract dark signals, yielding an image which is composed of primary x-rays and scatter 

20 x-rays; 

(g) producing, from said image D^, a low-resolution image Dfy representing said selected front detection locations; 

(h) acquiring a raw low-resolution image l^ from said selected rear detection locations and processing said image 
Irij to normalize it and to subtract dark signals, yielding an image D^, composed of substantially only said primary 
x-rays; 

25 (i) calculating a pair of low-resolution primary x-ray images D^p, and D^p, at said front detector assembly from a 

low-resolution dual-energy primary x-ray imaging pair composed of said image D,^, and said image D^; 
(j) calculating a low- resolution scatter x-ray image D fLSl at said front detector assembly by subtracting said image 
D^pi from said image and calculating a low-resolution scatter x-ray image D^, at said front detector assembly 
by subtracting said image D^p, from said image D^; 

30 (k) calculating a high -resolution scatter x-ray image D 1] j Sh by extending said low-resolution scatter x-ray image 

Djlsi t0 al1 of sa ' d front detection locations through interpolation and calculating a high-resolution scatter x-ray 
image D fHQh by extending said low-resolution scatter x-ray image DfHs, to all of said front detection locations through 
interpolation; and 

(I) calculating a high-resolution primary x-ray image D^ph at said front detector assembly by subtracting said image 
35 from said image D mh and calculating a high-resolution primary x-ray image DfLp h at said front detector 

assembly by subtracting said image D fLSh from said image Dfl_ h ; 

(m) whereby said images D^pp, and Dfl_p h form a high -resolution, two-dimensional, dual-energy primary x-ray 
image pair DfHp hJ of said subject after said scatter x-rays have been substantially eliminated, said image 
pair having a resolution substantially equal to the highest resolution available from said front detector assembly. 

40 

According to a fifth aspect of the present invention, there is provided a method for taking a two-dimensional x-ray 
image of a subject using a two-dimensional x-ray imaging system, said subject being composed substantially of two 
materials that interact differently with x-rays, said imaging system including, in physical sequence from front to back, 
an x-ray source, a front two-dimensional x-ray detector assembly having a plurality of detection locations, beam se- 

45 lection means, and a rear detector assembly assembly, said subject being between said x-ray source and said front 
detector assembly, said x-ray source being adapted to emit x-rays for passage through said subject, said x-rays in- 
cluding primary x-rays having their direction of travel unaltered by interaction with said subject and said x-rays including 
scatter x-rays having their direction of travel altered by interaction with said subject, said rear detector assembly in- 
cluding a rear low-energy two-dimensional x-ray detector having a plurality of selected rear low-energy detection lo- 

50 cations, a rear x-ray energy spectral filter, and a rear high-energy two-dimensional x-ray detector having a number of 
selected rear high-energy detection locations substantially equal to said plurality of selected rear low-energy detection 
locations, said beam selection means blocking the passage of substantially all of said scatter x-rays to said selected 
rear low-energy detection locations and said selected rear high -energy detection locations and permitting the passage 
of said primary x-rays to said selected rear low-energy detection locations and said selected rear high-energy detection 

55 locations, said selected rear low-energy detection locations and said selected rear high-energy detection locations 
receiving only primary x-rays, and selected front detection locations being those of said front detection locations inter- 
sected by x-ray projection lines extending from said x-ray source to said selected rear low-energy detection locations 
and said selected rear high-energy detection locations, said method comprising: 
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(a) illuminating said subject with said x-rays; 

(b) acquiring a raw high-resolution image l^ from said front detection locations and processing said image l fe to 
normalize it and to subtract dark signals, yielding a high-resolution image D m , which is composed of said primary 
x-rays and said scatter x-rays; 

s (c) producing, from said image D^, a low-resolution image D fl representing said selected front detection locations; 

(d) acquiring a raw low-resolution image l^, from said selected rear high-energy detection locations and processing 
said image I^j to normalize it and to subtract dark signals, yielding a low-resolution image D^i composed of 
substantially only said primary x-rays; 

(e) acquiring a raw low-resolution image I ^ from said selected rear low-energy detection locations and processing 
io said image 1^ to normalize it and to subtract dark signals, yielding a low-resolution image D^j composed of sub- 
stantially only said primary x-rays; 

(f) calculating a low-resolution primary x-ray image D m at said front detector assembly from a low-resolution dual- 
energy primary x-ray imaging pair composed of said image D^, and said image D^,; 

(g) calculating a low-resolution scatter x-ray image D^, at said front detector assembly by subtracting said low- 
is resolution primary x-ray image D fPJ from said low-resolution image D fl ; 

(h) calculating a high-resolution scatter x-ray image by extending said low-resolution scatter x-ray image D^, 
to all of said front detection locations through interpolation; and 

(i) calculating a high-resolution primary x-ray image Dfp h by subtracting said high-resolution scatter x-ray image 
Df£ h from said high -resolution image D^; 

20 (j) whereby said high-resolution primary x-ray image D fPh is a two-dimensional image of said subject after said 

scatter x-rays have been substantially eliminated, said image Dfp h having a resolution substantially equal to the 
highest resolution available from said front detector assembly. 

According to a sixth aspect of the present invention, there is provided a method for taking a two-dimensional x- 
25 ray image of a subject using a two-dimensional x-ray imaging system, said subject being composed substantially of 
two materials that interact differently with x-rays, said imaging system including, in physical sequence from front to 
back, an x-ray source, a front detector assembly, beam selection means, and a rear detector assembly, said subject 
being between said x-ray source and said front detector assembly, said x-ray source being adapted to emit x-rays for 
passage through said subject, said x-rays including primary x-rays having their direction of travel unaltered by inter- 
so action with said subject and said x-rays including scatter x-rays having their direction of travel altered by interaction 
with said subject, said front detector assembly including a front low-energy two-dimensional x-ray detector having a 
plurality of front low-energy detection locations, a front x-ray energy spectral filter, and a front high-energy two-dimen- 
sional x-ray detector having a number of front high-energy detection locations equal to said plurality of front low-energy 
detection locations, said rear detector assembly including a rear low-energy two-dimensional x-ray detector having a 
35 plurality of selected rear low-energy detection locations, a rear x-ray energy spectral filter, and a rear high-energy two- 
dimensional x-ray detector having a number of selected rear high-energy detection locations equal to said plurality of 
selected rear low-energy detection locations, said beam selection means blocking the passage of substantially all of 
said scatter x-rays to said selected rear low-energy detection locations and said selected rear high-energy detection 
locations and permitting the passage of said primary x-rays to said selected rear low-energy detection locations and 
40 said selected rear high -energy detection locations, said selected rear low-energy detection locations and said selected 
rear high-energy detection locations receiving only primary x-rays, selected front high-energy detection locations being 
those of said front high-energy detection locations intersected by x-ray projection lines extending from said x-ray source 
to said selected rear low-energy detection locations and said selected rear high-energy detection locations, selected 
front low-energy detection locations being those of said front low-energy detection locations intersected by said x-ray 
45 projection lines, said method comprising: 

(a) illuminating said subject with said x-rays; 

(b) acquiring a raw high-resolution image l^ from said front low-energy detection locations and processing said 
image l^ to normalize it and to subtract dark signals, yielding a high -resolution image D^, which is composed 

50 of said primary x-rays and said scatter x-rays; 

(c) producing, from said image Dflj,, a low-resolution image representing said selected front low-energy de- 
tection locations; 

(d) acquiring a raw high-resolution image l^ from said front high-energy detection locations and processing said 
image l Wh to normalize it and to subtract dark signals, yielding a high -resolution image D^, which is composed 

ss of said primary x-rays and said scatter x-rays; 

(e) producing, from said image D^, a low-resolution image D^, representing said selected front high-energy 
detection locations; 

(f) acquiring a raw low-resolution image l^, from said selected rear high-energy detection locations and processing 
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said image l^j to normalize it and to subtract dark signals, yielding a low-resolution image D^i composed of 
substantially only said primary x-rays; 

(g) acquiring a raw low-resolution image l^j from said selected rear low-energy detection locations and processing 
said image l^, to normalize it and to subtract dark signals, yielding a low-resolution image composed of sub- 

s stantially only said primary x-rays; 

(h) calculating a pair of low-resolution primary x-ray images D^p,; at said front high-energy detector and D^p, at 
said front low-energy detector from a low-resolution dual-energy primary x-ray imaging pair composed of said 
image D^, and said image D rU ; 

(i) calculating a low-resolution scatter x-ray image D^g, at said front high -energy detector by subtracting said 
10 image D^p, from said image D^, and calculating a low-resolution scatter x-ray image D^, at said front low-energy 

detector by subtracting said image D^p, from said image Dfy; 

(j) calculating a high-resolution scatter x-ray image by extending said low-resolution scatter x-ray image 
DfHsi to all of said front high-energy detection locations through interpolation and calculating a high-resolution 
scatter x-ray image Dfl_ Sh by extending said low-resolution scatter x-ray image Dfl_si to all of said front low-energy 

is detection locations through interpolation; and 

(k) calculating a high-resolution primary x-ray image D fHPh by subtracting said image DfH Sh from said image 
and calculating a high-resolution primary x-ray image Dfl_ Ph by subtracting said image Df,^ from said image D fLh ; 
(I) whereby said images D fHPh and Dfl_ Ph form a high-resolution, two-dimensional, dual-energy primary x-ray image 
pair DfHph, Dfl_ Ph of said subject after said scatter x-rays have been substantially eliminated, said image pair having 

20 a resolution substantially equal to the highest resolution available from said front detector assembly. 

The dual-energy x-ray data decomposition method in this invention consists in directly solving the dual energy 
nonlinear equation system, it includes: (1) Constructing an explicit quantitative equation system D^, = D H (b,s) and D L 
= D L (b,s) for each detector according to the nonlinear dual-energy x-ray imaging fundamental equation system in its 
25 original form. (2) reconstructing an equation system b = b(D H ,D L ) and s = s(D H ,D,J by numerically inverting the equation 
system of step 1.(3) Determining the desired values for b and s from the available data pair (Dh.DJ by using the 
numerical equations of step 2, or determining the desired values for D H ,D L , from the available data pair (b,s) by using 
the numerical equations of step 1 . (4) Maintaining the accuracy at each step to be as high as real number analytical 
solutions can provide. 

30 Embodiments of the present invention advantageously comprise an x-ray detection system that has the following 

features: (1) using a beam selection means to physically separates the primary x-rays from the scatter x-rays at a 
number of selected locations; (2) using a front detector assembly detects the primary and scatter x-rays, using a rear 
detector assembly detects only the primary x-rays at selected locations without scatter. 

The method for eliminating scatter x-rays advantageously includes: (1) using a dual-energy method to determine 
35 the primary image of the front detector at the corresponding selected locations from the image data of the rear detector 
assembly; (2) using an interpolation method to extend the scatter component of the front detector image determined 
at a the selected locations to the entire image area of the front detector. Following this, a complete separation of scatter 
and primary images for the front detector at high spatial resolution is achieved. 

Thus, embodiments of the present invention seek to provide an apparatus and method for substantially eliminating 
40 the effects of scatter on two-dimensional x-ray detectors. 

Embodiments of the present invention may also seek to provide an apparatus and a method for scatter-free dual- 
energy x-ray imaging using two-dimensional detectors and to provide two high accuracy material composition images 
of a subject at the spatial resolution as high as a two-dimensional detector array can provide. 

For a better understanding of the present invention, and to show how it may be carried into effect, reference shall 
45 now be made by way of example to the accompany drawings, in which: 

Fig. 1 is a diagram of the basic hardware of an embodiment of the present invention; 

Fig. 2 is a curve describing a typical x-ray source energy spectrum used in embodiments of the present invention; 
Fig. 3 is a flow diagram of the basic dual-energy data decomposition method using the hardware of Fig. 1 ; 
50 Fig. 4 is a flow diagram of the basic method to eliminate scatter using the hardware of Fig. 1 ; 

Fig. 5 is a diagram of a first embodiment of the present invention; 

Fig. 6 is a flow diagram of the method of the first embodiment using the hardware of Fig. 5; 
Figs. 7a and 7b are diagrams of a special case of the first embodiment of the present invention; 
Fig. 8 is a flow diagram of the method of the special case using the hardware of Figs. 7a and 7b; 
55 Fig. 9 is a diagram of a second embodiment of the present invention; 

Fig. 10 is a flow diagram of the method of the second embodiment using the hardware of Fig. 9; 
Fig. 11 is a diagram of a third embodiment of the present invention; 

Fig. 12 is a flow diagram of the method of the third embodiment using the hardware of Fig. 11; and 
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Figs. 13a to 1 3d is a graphically representation of a method for inverting the nonlinear dual-energy equation sys- 
tems. 

Introduction 

5 

The basic apparatus includes five components: (1) an x-ray source, (2) a front two-dimensional x-ray detector 
assembly, (3) a beam selector, (4) a rear two-dimensional x-ray detector assembly, and (5) a computer. 

As shown in Fig. 1 , the subject under examination 12 is located between the x-ray source 14 and the front detector 
assembly 16. The x-ray source 14 emits x-rays 30, preferably with an energy in the range of from 10 keV to 500 keV, 
10 as shown in Fig. 2. A portion of the x-rays 32 passes through the subject 1 2 directly to the front detector assembly 16 
without a change in their direction of propagation. These x-rays 32 are called the primary x-rays and convey true 
information about the subject 12. The remainder of the x-rays 34 are randomly scattered as a result of interaction with 
the material of the subject 1 2. These x-rays 34 are called scatter and cause a distortion of the true information. 

The front detector in the front detector assembly 1 6 has a number of detector types. The first uses thin film amor- 
's phous silicon as photodetection medium. An scintillation screen, is placed in close contact with the entire photosensitive 
area of the photodetector array. Preferably, the detector array has dimensions of 20 centimeters (cm) by 20 cm or 40 
cm by 40 cm for a single detector module. The cell size for this detector array is in the range of from approximately 60 
|im by 50 jim to approximately 1 mm by 1 mm. 

A second type of detector array uses an amorphous selenium film or selenium alloy film as the x-ray sensitive 
20 medium. Other typical two-dimensional detector arrays include but not limited to, charge-couple device (CCD) detec- 
tors, thin-film thallium-bromide-based detector arrays, avalanche silicon detector arrays, and phosphorstimulatable 
computed radiography screens. 

The front detector assembly 16 conveys an image that is a combination of the primary x-rays 32 and the scatter 
34, and is denoted by 

25 

D fh( x .y) = D fPh( x -y) + D fsh( x .y) m 

i 

where D f denotes an image in the front detector assembly 16 and (x,y) denotes the two-dimensional Cartesian 

30 coordinates of a cell of the front detector assembly 16. D fPh (x,y) denotes the contribution from the primary x-rays 32 
and D fSh (x,y) denotes the contribution from the scatter 34. 

The beam selector is sandwiched between the front detector assembly 16 and the rear detector assembly 26. A 
preferred embodiment of the x-ray beam selector 18 is a quantity of x-ray-absorbent material having a large number 
of holes 20. The holes 20 are fabricated such that their axes are aligned with the travel direction of the primary x-rays 

35 32. As a result, the holes 20 permit all x-rays traveling along the axes of the holes 20 to pass through, while all x-rays 
traveling in directions deviating slightly from the hole axes are completely absorbed by the bulk material of the beam 
selector 18. Because the holes 20 always have a finite size, a small portion of randomly scattered x-rays from the 
image subject 12 still can reach the rear detector assembly 26. However, as long as the hole size 20 is small and the 
thickness of the beam selector 18 is sufficiently large, this portion of the scatter 34 can be controlled to be negligibly 

40 small. Preferably, the holes 20 have a diameter that is in the range of approximately from 0.5 mm to 1 0 mm. Preferably, 
there are as many holes as possible in the beam selector 1 8. However, the beam selector material must occupy enough 
area so that the amount of scatter 34 is reduced to an insignificant quantity. A compromise based on these factors 
results in a pitch that is preferably between 2 mm and 50 mm. Preferably, the x-ray source 14 is located between 20 
cm and 150 cm from the rear surface 24 of the beam selector 18. The invention holds equally true when the x-ray 

45 source is a point source or has a finite size. 

The term "selected location' is defined as the location on the rear detector assembly 28, or on a rear detector, 
where, due to the function of the beam selector, only primary x-rays are received, and the scatter x-rays are substantially 
blocked. The "selected projection line" is defined as a straight line connecting the x-ray source to a point in the "se lected 
location'. The rear detector cells at the selected locations have a fixed geometric relation with some of the front detector 

50 cells. This relation is established by drawing a selected projection line from the x-ray source 14 through the beam 
selector 18 to the selected location. This selected projection line intersects the rear detector surface at a rear detector 
cell at a coordinate (i,j), and intersects the front detector surface at a front detector cell at a coordinate (x(i).y(j)). Here 
(x(i).yG)) denote the Cartesian coordinate (x.y) of the front detector cell in the front detector assembly 1 6 closest to the 
selected projection line. The relationship between (i,j) and (x(i),y(j)) is experimentally established for all of the holes 

55 20 of the beam selector 1 8 and stored in computer. The images represented by the composite of the signals from the 
detector cells of only one detector on the selected projection lines are low-resolution images and are represented by 
the subscript lower-case 1. The images represented by the composite of the signals from all the front detector cells 
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are high-resolution images and are represented by the subscript lower-case h. 

In connection with the material composition of the image subject, four quantities are defined. b(i,j) and s(i,j) are 
defined as the selected projection mass densities along the selected projection line (i,j). b(x,y) and s(x,y) are defined 
as the projection mass densities along the projection line (x,y). The "projection mass density" is defined as the integrated 
s total mass of the image subject along the projection line per unit area. Because the projection mass density is not 
dependent on the size of detector cells, b(x(i),y(j)) = b(i,j) and s(x(i),y(j)) = s(i,j). 

The method for eliminating scatter, shown in the flow diagram of Fig. 4, consists of the following steps: (1 ) illumi- 
nating the subject with x-rays from the x-ray source 14; (2) acquiring a low-resolution image pair l,Hi(ij) and l^flj) from 
the rear detector assembly and processing it to normalize and to subtract dark signals, yielding a low-resolution image 
10 pair Drf-i^i, j) and D^pj) that are functions of the subject materials, where "acquiring an image" is defined as transferring, 
via electronic control circuits, the electrical signals induced by the x-ray illumination on each detector cell from a detector 
array to computer memory; (3) solving the image pair D^O.j) and D rU (i,j) to determine the selected projection mass 
density of the subject materials b(i,j) and s(i,j) according to dual energy data decomposition methods; (4) acquiring a 
high-resolution image ^(x.yjfrom the front detector and processing it to normalize and to subtract dark signals, yielding 
is a high -resolution image Dfofx.y), which is the sum of primary x-rays and scatter x-rays; (5) calculating the low-resolution 
primary image of the front detector Dfp^xfO.yfl)) at the detector cells on the selected projection lines from the mass 
projection densities b(i,j) and s(i,j); (6) subtracting Djp^xpJ.y (j)) from D^fx^.yQ')) to calculate the low-resolution scatter 
component D fS ,(x(t) ( yQ)) of the image D m (x(i),yQ)) at the detector cells on selected projection lines; (7) interpolating 
D fSI (x(i)(x{j)) for those front detector cells not on the selected projection lines, yielding the high-resolution scatter image 
20 D^fx.y); and (8) subtracting the image D^ h (x,y) from D^x.y) to yield an image Dfp h (x,y), which is a full two-dimensional 
image of the subject at the front detector after scatter x-rays have been substantially eliminated. 

Because the data decomposition method of the present invention is rigorous and accurate and the mathematical 
calculations can be performed with good precision, as long as the experimental image data acquired from the detectors 
16, 26 is highly accurate, the final result will be a highly accurate primary image of the subject 12. 

25 

First embodiment 

In the first embodiment of the apparatus, shown in Fig. 5, the rear detector assembly 26 is constructed as a dual- 
energy x-ray imaging detector assembly. It has a low-energy two-dimensional detector 40, an x-ray energy spectral 

30 filter 42, and a high-energy two-dimensional detector 44. After the filter 42, x-ray energy is higher than that before the 
filter 42. The image detected by the low-energy detector 40 is denoted by D^iJ) and the image detected by the high- 
energy detector 44 is denoted by D^p.j). Preferably, the low-energy x-rays have an average energy of from 10 keV to 
1 00 keV and high-energy x-rays have an average energy of from 30 keVto 500 keV, with the high-energy x-rays having 
a higher energy than the low-energy x-rays. 

35 The purpose of measuring a pair of dual-energy primary images D^i.j) and D^i.j) of the rear detector assembly 

is to provide data for uniquely determining the corresponding primary low resolution image Dfp|(x(i),yG)) at the front 
detector. 

A flow diagram describing the method for determining a scatter-free image using the hardware of the first embod- 
iment is shown in Fig. 6. The x-ray source 14 emits x-rays with a uniform angular distribution and with an energy 
40 spectrum of O 0 (E). The x-rays passing through the subject 1 2 carry information on the thickness and material compo- 
sition of the subject 1 2, expressed as projection mass density in units of grams/centimeters 2 (g/cm 2 ). The image induced 
by the x-rays incident on the front detector 1 6 is denoted as D^x.y) and is 

45 D fh (x,y) = J [4> 0 (EJ x exp(-(ji b (E) x b(x,y) + (2) 

JLt s (E) x s(x,y))] x S*(E)dH + 
J <D S (E) x S f (E)dE 

50 

where u. b (E) is the mass absorption coefficient of bone tissue and |i s (E) is the mass absorption coefficient of soft 
tissue, with both n b (E) and u^(E) expressed in units of centimeter 2 /gram (cm 2 /g). Both of these values are known, 
having been determined experimentally and tabulated many years ago. Sf(E) is the x-ray spectral sensitivity (the elec- 
trical signal amplitude from the detector as a function of the number of x-rays with energy E incident upon the detector) 
55 of the front detector 16. The term J <D S (E) x Sf(E)dE represents the signal caused by scatter. 

The rear detector assembly 26 has two detectors 40, 44, so there are two low-resolution images D^i, j) and 
(i,j), which are 
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DrLi(i.j) = J [<S> 0 (E) x exp(-(ji b (E) x b(i,j) + (3a) 
m(E) x s (i, j ) ) ] x S rL (E)dE 

and 


E>rHi(i,j) = J [0> 0 (E) x exp(-(^i b (E) x b(ij) + (3b) 
p. s (E) X s(ij))] x S rH (E)dE 


Equations 3a and 3b constitute a simultaneous equation system, where the values for the signal pair D rlJ (i,j), D^, 
(i,j) are known quantities, being evaluated from the electrical signals of the rear detectors 40, 44. b(i,j) and s(i,j) are 
the unknown quantities for which equation pair 3a, 3b must be solved, as described below. Generally speaking, in 
mathematics, such a nonlinear simultaneous equation system may have an infinite number of solutions, may have 
multiple-value solutions, or may not have any solution. However, for the specific physical case of x-ray imaging, where 
the energy range is limited, preferably between 10 keV and 500 keV, it can be mathematically proved that a unique 
solution always exists for the nonlinear equation system 3a and 3b. 

Accurate b(i,j) and s(i,j) are calculated by the data decomposition method of the present invention, as will be 
described below. Now that the values for b(i,j) and s(i,j) are known, the front low- resolution scatter-free image Dfp,(x, 
y) can be obtained for those front detector cells (x(i).y(j)) that are on the selected projection lines. D fP ,(x(i),y(j)) is the 
signal induced by only primary x-rays at the detector cell (x(i),y(j)) on the front detector, and is 

D f? i(x(i) ,y(j) ) = I [O 0 (E) x exp(-{^i b (E) x (4) 

b(i, j) + M-s(E) x s(i, j ) ) ] x 
S f (E)dE 

where (x(i),y(j)) is the coordinate of the front detector cell (x,y) lying on the same selected projection line as the rear 
detector cell (i,j) and the energy dependent function <D 0 (E) x S^E) is given in calibration, as will be described in the 
data decomposition section below. 

Next, the low-resolution front scatter image D fS ,(x(i),y(j)) iS determined by applying equation 3, 

D«MD. (j)) = DflMO.yO) - D fp| (x(i),y(j)) 

Because of the scatter is largely caused by Compton scattering, which has a rather smooth angular distribution, 
the low-resolution scatter image D^fxfO.yQ)) can be extended to the entire (x.y) plane through interpolation without 
losing accuracy, yielding the high-resolution scatter image D^x.y). The high-resolution scatter image D^tx.y) is 
subtracted from the experimentally measured image D^x.y), yielding the high-resolution scatter-free signal Djp^x.y). 

It is now clear why it is crucial to acquire a pair of dual-energy x-ray images D^fiJ), D rH (i,j) to determine the scatter. 
The primary signals at the front detector D m (x(\),y{\)) can be uniquely determined only when the material composition 
image pair b(i,j) and s(i,j) is known (see equation 7). If only one image at the rear detector is used, it would always be 
found that the signal ratio between the front detector and the rear detector is dependent on the energy spectrum or 
dependent on the image subject. Consequently, the primary image on the front detector could not be determined. 

It is also noted that, as a supplement to this general rule, there are two special cases where the relationship 
between the signals on the image plane (i,j) and the image plane (x(i).y(j)) can be degenerated or simplified to be 
proportional to each other and independent of the image subject. Thus, for correcting scatter, it is not necessary to use 
the dual-energy method for these two special cases. 

The first special case occurs when the range of the x-ray energy spectral distribution is sufficiently narrow so that 
the x-rays can be approximated as having only a single energy or a wellndefined average energy E 0 . The signals on 
the front detector and the signals on the rear detector become a constant .independent of image subject and can be 
predetermined before using the system for imaging operations. 


9 


EP 0 782 375 A1 

The second special case is shown in Figs. 7a and 7b, where only one detector is used for acquiring both the high 
resolution image D^fx.y) and the low resolution image D ri (i,j). A mechanical device moves the beam selector to and 
away from the detectors front surface. When the beam selector is present, as in Fig. 7a, blocking the scatter from 
reaching the detector, a low resolution image D ri (i,j) is acquired. When the beam selector is absent, as in Fig. 7b, 

5 allowing all x-rays to reach the detector, a high resolution image D fh (x,y) is acquired, a portion of which is D fh (x(i),y(j). 
Here the subscript f denotes that the beam selector 18 is in its transmissive position not in front of the detector 16 and 
the subscript r denotes that the beam selector 18 is in its blocking position in front of the detector 16. In this special 
case, when the rest conditions are maintained substantially identical, if, and only if, the x-ray beam transmission through 
the holes is 100%, an energy-independent constant exists, which is independent of the image subject 12 and is pre- 

io determined when the image subject 12 is absent 

D fP1 (x(i),y(j)) = C g (i,j) x D ri (iJ) (5) 

15 Fig. 8 shows the data processing procedures for the second special case. 

Second embodiment 


20 


25 


35 


40 


45 


50 


In this embodiment, shown in Fig. 9, the front detector assembly 16 and the rear detector assembly 26 each have 
only one detector. The x-ray source 14 emits two consecutive pulses, a low-energy pulse followed by a high-energy 
pulse. 

As shown in Fig. 10, the acquired first image pair includes the high-resolution images from the front detector 16 
and are 


D eH h(x,y) = 1 [*oh(E) x exp(-(|i b (E) x b(x,y) + (6a) 
Ji s (E) x s(x,y))] x S £ (E)dE + 
30 ! 4> fs (E,x,y) x S e (E)dE 


and 


Df L h(x,y) = I [<E>ol(E) x exp(- (}i b .(E) x b(x,y) + (6b) 
fi s (E) x s (x,y) ) ] x S f (E)dE + 
/ O fs (E,x,y) x S f (E)dE 

and the second image pair includes the low-resolution images from the rear detector 26 and are 

D rK1 (i,j) = J [O 0H (E) x exp(-(ji b (E) x b(i,j) + (7a) 
|± s (E) x s (i, j ) ) ] x S r (E)dE 

and 


DrLi(i,j) = I [O 0L (E) x exp(-(u b (E) x b(i,j) + (7b) 
55 \i s (E) x s(i,j))] x S r (E)dE 

In the equation pair 7a, 7b, the acquired low-resolution image data are free of scatter radiation. By using the dual- 
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energy data decomposition methods described below, the simultaneous equation pair 7a, 7b is solved to find the 
solutions for the image pair of material composition b(i,j) and s(i,j). 

As described above with reference to the first embodiment, because the rear detector cell (i,j) and front detector 
cell (x(i).yQ)) lie on the same selected projection line, the low-resolution front detector primary image pair DfHp,(x(i),y 
(j)). Dfl_p,(x(i),yG)) can be determined from the rear detector primary image pair D^i,]), D^fiJ). The front detector 
scatter image pair D^M^yO)), D fl _ s ,(x(i) I y(j)) are found by the equations 

D|HBi«0.yO)) = D««l). y(i)) - D^MiJ.yQ)) (8a) 

and 

D tLSI (x(i),ya)) = D^fxfl.yU)) ■ D fLP| (x(i),yG)) (8b) 

As above, the low-resolution scatter images can be extended to all front detector cells not on the selected projection 
lines through interpolation without loss of accuracy to yield the high-resolution scatter image pair D fHSh (x,y), Dfl_ Sh (x, 
y). The high-resolution scatter-free images on the front detector assembly are denoted as D fHPh (x,y) and D^p^x.y) 
and are 

D iHPh( x ' y) = D fHh( x .y) - D fHSh( x -y) (9a) 

and 

D n.Ph( x -y) = Dfuifry) - D fi_sh< x .y) m 

The image pair D^p^y), D fLPh (x,y) is a pair of dual-energy x-ray images without scatter. This image pair in turn 
relates to the material composition of the subject by the equations 


D cH ph(x,y) = 1 [<X>oh(E) x exp(-(|i b (E) x b(x,y) + (10a) 
35 M-s(E) x s(x,y))l x S £ (E)dE 


and 


D fL Ph(x,y) = 1 [*ol(E) x exp(-(u b (E) x b(x,y) + (10b) 
|1 S (E) x s (x,y) ) ] x S f (E)dE 


Thus, in addition to providing one scatter-free image, this embodiment provides a pair of scatter-free dual-energy 
images in the equation pair 10a, 10b. This equation pair is the fundamental dual-energy x-ray imaging equation system 
with the unprecedented feature that scatter radiation has been essentially removed from the two-dimensional detector. 
In the equation pair 10a, 10b, the values of O^p^y) and D (HPh (x,y) are known from the above described calculations 
conducted on the image pair D^x.y) D^fx.y) directly measured from the front detectors 16, and on the image pair 
D rij(U). CWU) directly measured from the rear detector 26. The unknown values are the two material composition 
images b(x,y) and s(x,y). 

The dual-energy x-ray data decomposition method can be further applied to the equation pair 1 0a, 1 0b. As a result, 
by using the quantitative relationships b = b(D H ,D L ) and s = s(D H ,D L ) provided by the data decomposition method, a 
pair of high-resolution images b(x,y) and s(x,y) are readily obtained point by point for all front detector cells (x,y). The 
solution of the two-component material composition images b(x.y) and s(x,y) has a spatial resolution as high as the 
front detector 16 can provide. 

An alternate to the second embodiment substitutes an x-ray source generates x-rays continuously, alternating 
between high-energy x-rays and low-energy x-rays. 
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. Another alternate to the second embodiment inserts an x-ray energy filter between the x-ray source and the subject 
at the moment the x-ray source switches to generate the high-energy x-rays. The synchronization between the insertion 
of the filter and the switching high-voltage or double-pulse is preferably implemented by using a motor drive. The filter 
absorbs more of the low-energy x-rays, resulting in an increase in the energy difference between the low-energy x- 
5 rays and the high-energy x-rays. 

This invention not only provides a method and apparatus for removing scattering from two-dimensional detectors, 
but at the same time also provides a method and apparatus for scatter-free dual-energy x-ray imaging using two- 
dimensional detectors. 

The interrelationship between the method for removing scatter radiation in two-dimensional detectors and the 
10 method for dual-energy x-ray imaging using two-dimensional detectors can be summarized as follows: 

1. The method for removing scatter radiation from two-dimensional detectors utilizes and hinges on the method 
of dual-energy x-ray imaging free of scatter. Without dual-energy x-ray imaging, the scatter radiation cannot be 
accurately removed. 

15 2. The method for dual-energy x-ray imaging using two-dimensional detectors utilizes and hinges on the method 

of removing scatter radiation. Without substantially removing scatter from two-dimensional detectors, the accuracy 
of dual-energy x-ray imaging would be so degraded as to be meaningless. This invention solves both problems in 
a unified system. 

20 Third Embodiment 

The third embodiment, uses a pair of two-dimensional detector assemblies. The x-ray source 14 is a constant 
potential x-ray source that emits steady state x-rays, single pulse x-rays, or repetitive pulse x-rays with the same energy 
spectrum. The front detector assembly 16 has a low-energy two-dimensional detector 50, an x-ray energy spectral 

25 filter 52, and a high-energy two-dimensional detector 54. The rear detector assembly 26 also has a low-energy two- 
dimensional detector 56, an x-ray energy spectral filter 58, and a high-energy two-dimensional detector 60. The filters 
52, 58 operate in the conventional manner as described above with reference to the first embodiment. The front high- 
energy detector 54 is sensitive to higher x-ray energies than the front low-energy detector 50 and the rear high-energy 
detector 60 is sensitive to higher x-ray energies than the rear low-energy detector 56. 

30 This embodiment requires only single-pulse, constant-energy x-rays, as in the first embodiment, rather than the 

dual-energy x-ray source of the second embodiment As shown in Fig. 12, following a single-pulse x-ray, two pairs of 
dual-energy x-ray images are acquired. The first pair includes the high-resolution images Df^x.y) from the low-energy 
front detector 50 and D^x.y) from the high-energy front detector 54. 

The second pair includes the low-resolution images D^O.j) from the low-energy rear detector 56 and D rH ,(i,j) from 

35 the high-energy rear detector 60. 

As described above with relation to the second embodiment, the image pair D rU (i,j) and D^fl.j) does not contain 
scatter. Thus, a pair of material composition images b(i,j), s(i,j) can be calculated for each cell in the (i,j) plane. This 
image pair is used to first determine the low-resolution primary image pair D fHP ,(x(i) l y(j)), D^p^xfO.yO)) and then to 
determine scatter image pair DjHs^O.yO)), DftjsrfxtiJ.yfl)) in the same way as for the equation pair 14a, 14b. Scatter 

40 image pair D^^x^.yO)), DfLsiMO.yG)) is used to further determine the high-resolution scatter image pair D^^y), 
D fLSh( x >y)> as described above in the second embodiment. By subtracting the calculated high-resolution scatter image 
P air D fHSh( x >y). D fLSh( x .y) from the front detector 16 high-resolution image pair D fHh (x,y), D^rx.y), the high spatial 
resolution image pair free of scatter D fHPh (x,y) and D^p^y) is numerically obtained. Using the inversion method 
described below, a pair of high -accuracy, high-resolution material composition images b(x,y) and s(x,y) can be further 

<5 obtained. 

Data Decomposition Method 

The following is a step-by-step description of the data decomposition method summarized above and shown in 
so Fig. 3. 

The first step is to construct the two simultaneous numerical equations D L = D L (b,s) and D H = D H (b,s) in three- 
dimensional space. A preferred method for doing this is to determine the function S(E) through the standard absorption 
method. 

An absorption curve is measured by using a collimated narrow primary x-ray beam. The electrical signal from a 
55 single detector cell D(t) as a function of the absorption plate thickness t is experimentally determined and is related to 
S(E) through the equation 
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D(t)= J S(E)X O 0 (E) X exp(-u.(E) X t)dE (1 ^ 

Since the mass absorption coefficient \i(B) of the absorption plate material is known, the function S(E) can be 
determined. This method is especially convenient for the internal conversion type of two-dimensional x-ray detectors. 
In these detectors, the detection efficiency and detector energy response function can be expressed in a simple ana- 
lytical expression with few unknown parameters to be solved, 

S(E) = {S 0 [l-exp(-(i 0 (E) x d) ] x aE} x (12 ) 
exp{-(l 1 (E) x di - fl 2 (E) x d 2 ) 

where S 0 (E)= [l-expf-u^E) X d)] x aE is the electrical signal amplitude induced by x-ray photons with energy E, 
(E) is the mass absorption coefficient of the detector's conversion layer, d is the thickness of the conversion layer of 
the detector cell, and where S^E) = exp(-u.-,(E) X 6, - u^(E) X d 2 ) is the x-ray transmission after leaving the image 
subject to the detector surface, m (E) and jj^E) are the absorption coefficients of two given materials, and and d 2 
are the thickness values of these materials. 

When the x-ray energy spectrum <X> 0 (E) is separately measured, these unknown parameters a, d, d 1f and d 2 are 
determined by using standard least square parameter-fitting techniques through equation 12. Then, the function S(E) 
is obtained to a high degree of accuracy for a single cell. 

Once the value for S(E) is determined to the desired accuracy, the dual-energy signals as a function of the material 
composition of the subject can be calculated through the equations 

D H = \ S(E) X <P 0 (E)x exp(-(jx b (E) x b + (13a) 
(I S (E) X s)dE 

and 


D L = | S(E) X 4> 0 (E)x exp(-(|i b (E) X b + (13b) 
|i s (E) X s)dE 


where u, b (E) and u, s (E) are the well-documented mass absorption coefficients for bone tissue and soft tissue, respec- 
tively. The mass surface densities b and s are assigned values that sufficiently cover the real range of the subject 12. 

Another preferred method for constructing the quantitative explicit functions D L = D L (b,s) and D H = D H (b,s) is to 
conduct direct measurements of signals D L and D H at a number of selected b and s values. The entire functions Dl = 
D L (b,s) and D H = D H (b,s) are obtained from the directly measured data points by using standard two-dimensional 
interpolation algorithms. The interpolation in this case is valid because the functions D L = D L (b,s) and D H = D H (b,s) 
are continuous, smooth, and monotonous. 

The second step is to determine the material composition images b and s as functions of the image pair D H , D L . 
The procedures for obtaining a simultaneous equation system for b(D H , D L ) and s(D H , D L ) are shown graphically in Figs. 
1 3a to 1 3d. To do so, the simultaneous equation pair 24a, 24b must be inverted. A preferred method of inversion is as 
follows: (1) as in Figs. 13a and 13b, assign a pair of values in the desired range to b and s corresponding to one of 
the coordinate points in the (b,s) plane so that b = b n , and s = s m , where n = 0,1,2.. .N, and m = 0,1,2.. .M. Typical N 
and M values are in the range of between approximately 50 and approximately 5,000. The larger N and M, the higher 
the accuracy of the results. However, the largest values for N and M are limited by the available capacity of computer 
memory and computing speed. From the two numerical equations representing the three-dimensional surfaces D H (b, 
s) and D L (b,s), determine a pair of D H and D L values so that D H [n,m] = D H (b = b n , s = s m ) and DJn.m] = D L (b = b n , s 
= sj, where D H [n,m] and DJn.m] are two specific real numbers, and (2) as in Figs. 13c and 13d, replot the four 
numbers D H [n,m], DJn.m], b n , and s m to provide a pair of data points on the three-dimensional surfaces b(D H ,D L ) and 
s ( d H' D l)- Tn9 data P oint on tne three-dimensional surface b(D H ,D L ) is D H = D H [n,m], D L = DJn.m], b = b n , and the 
data point on the three-dimensional surface s(D H ,D L ) is D H = D H [n,m], D L = DJn.m], s = s m . After going through all the 
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b = b n values (b 0 ,b v b 2> ...,b N ) and all the s = s m values (S 0 ,S 1 ,S 2 ,...,S M ), the essential part of the inversion task is 
complete. However, for the purpose of storing the inverted arrays b = b(D H ,D L ) and s = S(D H ,D L ), the step sizes of D H 
= D H [n,m] and D L = DJn.m] must be adjusted. In the inverted space, D H and D L are basis coordinates. From the N x 
M data points, only J data points are selected for D H and only K data points are selected for D L , where J and K are 
approximately in the same range as N and M. In the final form after the second step, two two-dimensional arrays are 
obtained and stored: b = b(D H ,D L ) and s = s (D H ,D L ), where D H = D H [j], D L = DJk]; j=0,1 t 2,... f J, D H [j] > D H [j+1] and 
k=0 l 1,2 l ...,K, DJk] > DJk+1]. Two additional one-dimensional arrays D H [j] and DJk] are also stored. Arrays D H [j] and 
DJk] are saved so that accuracy as high as real number calculations can provide is maintained. 

The third step is to find the desired results from the input data according to the established equations. The desired 
values for b and s at each cell location is determined by inserting the available data pair (Dh.Dl) into the numerical 
equations of step 2. Conversely, the desired values for D H , D L , or only one of them if only one is needed, at each 
discrete cell location is determined by inserting the available data pair (b,s) into the numerical equations of step 1 . 

The final step is to maintain the accuracy of the values for b and s in order to maintain a continuous domain function. 
The following procedures ensure elimination of the errors in connection with finite steps in data processing. 

In step 1, in the process of constructing the equation pair for D H [n,m] = D H (b = b n , s = s m ) and DJn.m] = D L (b = 
b n» s = s m). for eacn P air of values of b n and s m , the D H [n,m] and DJn.m] are measured or calculated to an accuracy 
of real numbers. D H [n,m] and DJn.m] are stored in computer as real number arrays. 

In step 2, the inversion process, including replotting in D H space and D L space, introduces no errors due to the 
data processing. The step sizes can be changed without losing any accuracy as long as values for D^, = D H [j] are 
selected that are exactly equal to one of the D H [n,m] values that satisfies the condition D H [j-1] > D H [j] > D H 0+1], and 
values for D L = DJk] are selected that are exactly equal to one of the DJn.m] values that satisfies the condition D Jk- 
1]>DJk]>DJk+1]. 

In step 3, for each measured dual-energy signal data pair (D HEX> D LEX ), first to find out the closest j and k values 
according to the criteria: D H [j] > D HEX > D H [j+1 ] and DJk] > D^x > DJk+1 ). From the index values j and k, the closest 
b and s are first determined as b 0 = b 0 (D H [j], DJk]) and s 0 = s 0 (D H [j], DJk]). The following equations give b and s 
values to an accuracy as high as real number calculations can provide: 

b = b 0 (D H [ j] ,D L [k] ) + [3b(D H ,D L ) /3D L ] DK=DH (j]; dl=dujc] X (14a) 

(D LEX - D L [k]) + [ 3b { Dh / Dl ) / 3Dh ] DH=DH ( j ] ; DL=DL ( k] X 

(Dhex - DmIj]) + higher order terms 

and 

S = S 0 (D H [ j] ,D L [k] ) + [3s(D H ,D L )/8D L ] D H= D H[j]; DL=DL[)C) x (14b) 
(D LEX - D L [k]) + [3s (D H/ D L ) /3D H ] D H=DH{j]; DL=DL{k) X 

(Dhex - D H [ j ] ) + higher order terms 

where the values for the higher order terms are found in standard calculus textbooks. 

Also in step 3, if the image pair D H and D L from a given material composition data pair (b ex> s ex ) is to be found, D H 
and D L is obtained to an accuracy of real numbers by using similar standard Taylor expressions. 

Thus, the procedures described above provide methods for directly solving the nonlinear dual-energy x-ray imaging 
fundamental equation systems in its original form with arbitrary bremsstrah lung spectra at an accuracy as high as using 
real number analytical calculations can provide. 

The following is a list of possible further variations in the embodiments: 

(1) Some well-established computation tools, such as sorting algorithms or database procedures, can be used to 
carry out the inversion process described above. The use of any such software package to conduct the above- 
described inversion process is contemplated by the present invention. 

(2) All of the steps described above, including the data decomposition method and the scatter elimination method, 
can be combined together to various degrees, from combining any two steps to combining all the steps into one 
procedure. For example, image pair (D^.D^) can be used directly to find (Dfl_|p,Dfl_p) without explicitly determining 
(b,s). One way of doing so is to construct a pair of quantitative relationships D^p = (D^.D^) and D^ = (D^.D^ 
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in a data base and storing them. Any such combination methods are contemplated by the present invention. 

It is intended that the scope of the invention be limited not by this detailed description, but rather by the claims 
appended hereto. 


Claims 

1 . A two-dimensional x-ray imaging system for taking images of a subject, said system comprising: 

(a) an x-ray source adapted to emit x-rays for passage through said subject, said x-rays including primary x- 
rays having their direction of travel unaltered by interaction with said subject and said x-rays including scatter 
x-rays having their direction of travel altered by interaction with said subject; 

(b) x-ray detection means for detecting said x-rays in two dimensions; 

(c) said subject being located between said x-ray source and said x-ray detection means; 

(d) said detection means including high-resolution detection means for detecting said primary x-rays and said 
scatter x-rays simultaneously; 

(e) beam selection means for blocking passage of substantially all of said scatter x-rays and permitting passage 
of a portion of said primary x-rays; 

(f) said detection means including tow-resolution detection means for detecting only said x-rays passed by 
said beam selection means. 

2. The x-ray imaging system of claim 1 wherein said beam selection means is substantially composed of an x-ray- 
absorbent material having a plurality of holes, the axes of said holes being parallel to the direction of travel of said 
primary x-rays. 

3. The x-ray imaging system of claim 1 wherein said detection means includes a detector assembly and wherein said 
beam selection means has two positions, a transmissive position wherein said high-resolution detection means 
includes said detector assembly, and a blocking position wherein said low-resolution detection means includes 
said detector assembly. 

4. The x-ray imaging system of claim 1 wherein said high-resolution detection means includes a front detector as- 
sembly and said low-resolution detection means includes a rear detector assembly, and wherein said system 
comprises, in physical sequence from front to rear, said x-ray source, said front detector assembly, said beam 
selection means, and said rear detector assembly, said subject being located between said x-ray source and said 
front detector assembly. 

5. The x-ray imaging system of claim 4 wherein said x-ray source is adapted to emit two x-ray pulses of different 
average energies. 

6. The x-ray imaging system of claim 4 wherein said rear detector assembly includes, in physical sequence from 
front to rear, a rear low-energy detector, a rear x-ray energy spectral filter, and a rear high-energy detector, said 
rear low-energy detector including a plural ity of x-ray-sensitive detector cells arranged in a substantially rectangular 
matrix, said rear high-energy detector including a plurality of x-ray-sensitive detector cells, the arrangement and 
quantity of said rear high-energy detector cells being substantially the same as the arrangement and quantity of 
said rear low-energy detector cells. 

7. The x-ray imaging system of claim 6 wherein said front detector assembly includes, in physical sequence from 
front to rear, a front low-energy detector, a front x-ray energy spectral filter, and a front high-energy detector, said 
front low-energy detector including a plurality of x-ray-sensitive detector cells arranged in a substantially rectan- 
gular matrix, said front high-energy detector including a plurality of x-ray-sensitive detector cells, the arrangement 
and quantity of said front high-energy detector cells being substantially the same as the arrangement and quantity 
of said front low-energy detector cells. 

8. A method for taking a two-dimensional x-ray image of a subject using a two-dimensional x-ray imaging system, 
said subject being composed substantially of two materials that interact differently with x-rays, said imaging system 
including an x-ray source, beam selection means, and a two-dimensional x-ray detector assembly having a plurality 
of detection locations, said subject being between said x-ray source and said detector assembly, said x-ray source 
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being adapted to emit x-rays for passage through said subject, said x-rays including primary x-rays having their 
direction of travel unaltered by interaction with said subject, said x-rays including scatter x-rays having their direc- 
tion of travel altered by interaction with said subject, said beam selection means having a blocking position wherein 
passage of substantially all of said scatter x-rays to said detector assembly is blocked and passage of said primary 
x-rays to said detector assembly is permitted, said beam selection means having a transmissive position wherein 
passage of substantially all of said primary x-rays and said scatter x-rays to said detector assembly is permitted, 
and selected detection locations being those of said detection locations receiving only said primary x-rays when 
said beam selection means is in said blocking position, said method comprising: 

(a) illuminating said subject with said x-rays; 

(b) moving said beam selection means to said transmissive position; 

(c) acquiring a raw high-resolution image l^ from said detection locations and processing said image l m to 
normalize it and to subtract dark signals, yielding a high -resolution image D m composed of said primary x- 
rays and said scatter x-rays; 

(d) producing, from said image D^, a low- resolution image D f , representing said selected detection locations; 

(e) moving said beam selection means to said blocking position; 

(f) acquiring a raw low-resolution image l ri from said selected detection locations and processing said image 
l r , to normalize it and to subtract dark signals, yielding a low-resolution image D ri composed of substantially 
only said primary x-rays; 

(g) calculating a low-resolution primary x-ray image Dfp, by multiplying said image D ri by a predetermined 
constant that is independent of said subject; 

(h) calculating a low-resolution scatter x-ray image D e , by subtracting said image D^, from said image D fl ; 

(i) calculating a high-resolution scatter x-ray image by extending said low-resolution scatter x-ray image 
Dfs, to all of said detection locations through interpolation; and 

(j) calculating a high-resolution primary x-ray image Dfp h by subtracting said high-resolution scatterx-ray image 
Dfs h from said high-resolution image D^; 

(k) whereby said image Dfp h is a high-resolution, two-dimensional primary x-ray image of said subject after 
said scatter x-rays have been substantially eliminated, said image Dfp h having a resolution substantially equal 
to the highest resolution available from said detector assembly. 

A method for taking a two-dimensional x-ray image of a subject using a two-dimensional x-ray imaging system, 
said subject being composed substantially of two materials that interact differently with x : rays, said imaging system 
including, in physical sequence from front to back, an x-ray source, a front two-dimensional x-ray detector assembly 
having a plurality of front detection locations, beam selection means, and a rear two-dimensional x-ray detector 
assembly having a plurality of selected rear detection locations, said subject being between said x-ray source and 
said front detector assembly, said x-ray source being adapted to emit x-rays for passage through said subject, 
said x-rays having an energy spectrum that can be approximated as having a single energy Eq, said x-rays including 
primary x-rays having their direction of travel unaltered by interaction with said subject and said x-rays including 
scatter x-rays having their direction of travel altered by interaction with said subject, said beam selection means 
blocking the passage of substantially all of said scatter x-rays to said selected rear detection locations and per- 
mitting the passage of said primary x-rays to said selected rear detection locations, said selected rear detection 
locations receiving only primary x-rays, and selected front detection locations being those of said front detection 
locations intersected by x-ray projection lines extending from said x-ray source to said selected rear detection 
locations, said method comprising: 

(a) illuminating said subject with said x-rays; 

(b) acquiring a raw high-resolution image l^ from said front detection locations and processing said image l m 
to normalize it and to subtract dark signals, yielding a high -resolution image D m composed of primary x-rays 
and scatter x-rays; 

(c) producing, from said image D m> a low-resolution image D fl representing said selected front detection loca- 
tions; 

(d) acquiring a raw low-resolution image l rt from said selected rear detection locations and processing said 
image l rt to normalize it and to subtract dark signals, yielding a low-resolution image D ri composed of substan- 
tially only said primary x-rays; 

(e) calculating a low-resolution primary x-ray image D^, at said front detector assembly by multiplying said 
image D r) by a predetermined constant that is independent of said subject; 

(f) calculating a low-resolution scatter x-ray image Dfs, at said front detector assembly by subtracting said 
image D fR from said image D fi ; 
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(g) calculating a high-resolution scatter x-ray image D,s h by extending said low-resolution scatter x-ray image 
Dfsj to all of said front detection locations through interpolation; and 

(h) calculating a high-resolution primary x-ray image Dfp h by subtracting said high-resolution scatter x-ray 
image Dfg h from said high-resolution image D m ; 

(i) whereby said image D fPh is a high-resolution, two-dimensional primary x-ray image of said subject after 
said scatter x-rays have been substantially eliminated, said image D,p h having a resolution substantially equal 
to the highest resolution available from said front detector assembly. 

10. A method for taking a two-dimensional x-ray image of a subject using a two-dimensional x-ray imaging system, 
said subject being composed substantially of two materials that interact differently with x-rays, said imaging system 
including, in physical sequence from front to back, an x-ray source, a front two-dimensional x-ray detector assembly 
having a plurality of front detection locations, beam selection means, and a rear two-dimensional x-ray detector 
assembly having a plurality of selected rear detection locations, said subject being between said x-ray source and 
said front detector assembly, said x-ray source being adapted to emit x-rays at two different average energy levels 
for passage through said subject, said x-rays including primary x-rays having their direction of travel unaltered by 
interaction with said subject and said x-rays including scatter x-rays having their direction of travel altered by 
interaction with said subject, said beam selection means blocking the passage of substantially all of said scatter 
x-rays to said selected rear detection locations and permitting the passage of said primary x-rays to said selected 
rear detection locations, said selected rear detection locations receiving only primary x-rays, and selected front 
detection locations being those of said front detection locations intersected by x-ray projection lines extending from 
said x-ray source to said selected rear detection locations, said method comprising: 

(a) illuminating said subject with x-rays of a first average energy level H; 

(b) acquiring a raw high-resolution image l^ from said front detection locations and processing said image 
IfHh t0 normalize it and to subtract dark signals, yielding an image D Mh which is composed of primary x-rays 
and scatter x-rays; 

(c) producing, from said image D mh , a low-resolution image representing said selected front detection 
locations; 

(d) acquiring a raw low-resolution image l^, from said selected rear detection locations and processing said 
image l^, to normalize it and to subtract dark signals, yielding an image D^, composed of substantially only 
said primary x-rays; 

(e) illuminating said subject with x-rays of a second average energy level L; 

(f) acquiring a raw high-resolution image from said front detection locations and processing said image 
Iflj, to normalize it and to subtract dark signals, yielding an image which is composed of primary x-rays 
and scatter x-rays; 

(g) producing, from said image Df^, a low-resolution image Dfy representing said selected front detection 
locations; 

(h) acquiring a raw low-resolution image l^, from said selected rear detection locations and processing said 
image l^, to normalize it and to subtract dark signals, yielding an image composed of substantially only 
said primary x-rays; 

(i) calculating a pair of low-resolution primary x-ray images D^p, and D^p, at said front detector assembly 
from a low-resolution dual-energy primary x-ray imaging pair composed of said image D^, and said image D^,; 
(j) calculating a low-resolution scatter x-ray image D^, at said front detector assembly by subtracting said 
image D^p, from said image D^, and calculating a low-resolution scatter x-ray image D^s, at said front detector 
assembly by subtracting said image D fHPj from said image D^; 

(k) calculating a high-resolution scatter x-ray image by extending said low-resolution scatter x-ray image 
D fLSi to a " of said f ront detection locations through interpolation and calculating a high-resolution scatter x-ray 
image D,^ by extending said low-resolution scatter x-ray image D^q, to all of said front detection locations 
through interpolation; and 

(I) calculating a high-resolution primary x-ray image D 1H p h at said front detector assembly by subtracting said 
image D msh from said image and calculating a high-resolution primary x-ray image Dfl_p h at said front 
detector assembly by subtracting said image from said image D^; 

(m) whereby said images D^ph and Dfl_ Ph form a high-resolution, two-dimensional, dual-energy primary x-ray 
image pair D^, Dfl_ Ph of said subject after said scatter x-rays have been substantially eliminated, said image 
pair having a resolution substantially equal to the highest resolution available from said front detector assembly 

11. A method for taking a two-dimensional x-ray image of a subject using a two-dimensional x-ray imaging system, 
said subject being composed substantially of two materials that interact differently with x-rays, said imaging system 
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including, in physical sequence from front to back, an x-ray source, a front two-dimensional x-ray detector assembly 
having a plurality of detection locations, beam selection means, and a rear detector assembly assembly, said 
subject being between said x-ray source and said front detector assembly, said x-ray source being adapted to emit 
x-rays for passage through said subject, said x-rays including primary x-rays having their direction of travel unal- 
tered by interaction with said subject and said x-rays including scatter x-rays having their direction of travel altered 
by interaction with said subject, said rear detector assembly including a rear low-energy two-dimensional x-ray 
detector having a plurality of selected rear low-energy detection locations, a rear x-ray energy spectral filter, and 
a rear high-energy two-dimensional x-ray detector having a number of selected rear high-energy detection loca- 
tions substantially equal to said plurality of selected rear low-energy detection locations, said beam selection means 
blocking the passage of substantially all of said scatter x-rays to said selected rear low-energy detection locations 
and said selected rear high-energy detection locations and permitting the passage of said primary x-rays to said 
selected rear low-energy detection locations and said selected rear high-energy detection locations, said selected 
rear low-energy detection locations and said selected rear high-energy detection locations receiving only primary 
x-rays, and selected front detection locations being those of said front detection locations intersected by x-ray 
projection lines extending from said x-ray source to said selected rear low-energy detection locations and said 
selected rear high-energy detection locations, said method comprising: 

(a) illuminating said subject with said x-rays; 

(b) acquiring a raw high-resolution image l fh from said front detection locations and processing said image l^ 
to normalize it and to subtract dark signals, yielding a high-resolution image D^, which is composed of said 
primary x-rays and said scatter x-rays; 

(c) producing, from said image D m , a low-resolution image D fl representing said selected front detection loca- 
tions; 

(d) acquiring a raw low-resolution image l rH , from said selected rear high-energy detection locations and 
processing said image l^, to normalize it and to subtract dark signals, yielding a low-resolution image D^i 
composed of substantially only said primary x-rays; 

(e) acquiring a raw low-resolution image l^, from said selected rear low-energy detection locations and 
processing said image l^ to normalize it and to subtract dark signals, yielding a low-resolution image D^, 
composed of substantially only said primary x-rays; 

(f) calculating a low-resolution primary x-ray image Dfp, at said front detector assembly from a low-resolution 
dual-energy primary x-ray imaging pair composed of said image D^i and said image D^; 

(g) calculating a low-resolution scatter x-ray image D,s, at said front detector assembly by subtracting said 
low-resolution primary x-ray image D^, from said low-resolution image D fl ; 

(h) calculating a high-resolution scatter x-ray image by extending said low-resolution scatter x-ray image 
Djg, to all of said front detection locations through interpolation; and 

(i) calculating a high-resolution primary x-ray image Dfp h by subtracting said high-resolution scatterx-ray image 

from said high-resolution image D^, 
(j) whereby said high-resolution primary x-ray image D fPh is a two-dimensional image of said subject after said 
scatter x-rays have been substantially eliminated, said image D fPh having a resolution substantially equal to 
the highest resolution available from said front detector assembly. 

12. A method for taking a two-dimensional x-ray image of a subject using a two-dimensional x-ray imaging system, 
said subject being composed substantially of two materials that interact differently with x-rays, said imaging system 
including, in physical sequence from front to back, an x-ray source, a front detector assembly, beam selection 
means, and a rear detector assembly, said subject being between said x-ray source and said front detector as- 
sembly, said x-ray source being adapted to emit x-rays for passage through said subject, said x-rays including 
primary x-rays having their direction of travel unaltered by interaction with said subject and said x-rays including 
scatter x-rays having their direction of travel altered by interaction with said subject, said front detector assembly 
including a front low-energy two-dimensional x-ray detector having a plurality of front low-energy detection loca- 
tions, a front x-ray energy spectral filter, and a front high-energy two-dimensional x-ray detector having a number 
of front high-energy detection locations equal to said plurality of front low-energy detection locations, said rear 
detector assembly including a rear low-energy two-dimensional x-ray detector having a plurality of selected rear 
low-energy detection locations, a rear x-ray energy spectral filter, and a rear high-energy two-dimensional x-ray 
detector having a number of selected rear high-energy detection locations equal to said plurality of selected rear 
low-energy detection locations, said beam selection means blocking the passage of substantially all of said scatter 
x-rays to said selected rear low-energy detection locations and said selected rear high-energy detection locations 
and permitting the passage of said primary x-rays to said selected rear low-energy detection locations and said 
selected rear high-energy detection locations, said selected rear low-energy detection locations and said selected 
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rear high-energy detection locations receiving only primary x-rays, selected front high-energy detection locations 
being those of said front high-energy detection locations intersected by x-ray projection lines extending from said 
x-ray source to said selected rear low-energy detection locations and said selected rear high-energy detection 
locations, selected front low-energy detection locations being those of said front low-energy detection locations 
intersected by said x-ray projection lines, said method comprising: 

(a) Illuminating said subject with said x-rays; 

(b) acquiring a raw high-resolution image l^ from said front low-energy detection locations and processing 
said image l^ to normalize it and to subtract dark signals, yielding a high-resolution image D 0 , which is 
composed of said primary x-rays and said scatter x-rays; 

(c) producing, from said image D^, a low-resolution image D fL , representing said selected front low-energy 
detection locations; 

(d) acquiring a raw high-resolution image l^ from said front high-energy detection locations and processing 
said image l^ to normalize it and to subtract dark signals, yielding a high-resolution image D^, which is 
composed of said primary x-rays and said scatter x-rays; 

(e) producing, from said image D^, a low-resolution image D^,, representing said selected front high-energy 
detection locations; 

(f) acquiring a raw low-resolution image l^, from said selected rear high-energy detection locations and 
processing said image l^, to normalize it and to subtract dark signals, yielding a low-resolution image D^, 
composed of substantially only said primary x-rays; 

(g) acquiring a raw low-resolution image l^, from said selected rear low-energy detection locations and 
processing said image l^, to normalize it and to subtract dark signals, yielding a low-resolution image D^, 
composed of substantially only said primary x-rays; 

(h) calculating a pair of low-resolution primary x-ray images D^p, at said front high-energy detector and D^p, 
at said front low-energy detector from a low-resolution dual-energy primary x-ray imaging pair composed of 
said image D^i and said image D^; 

(i) calculating a low-resolution scatter x-ray image D^s, at said front high-energy detector by subtracting said 
image D fHP , from said image and calculating a low-resolution scatter x-ray image D^, at said front low- 
energy detector by subtracting said image D^p, from said image D^; 

(j) calculating a high-resolution scatter x-ray image by extending said low-resolution scatter x-ray image 
Dfns, to all of said front high-energy detection locations through interpolation and calculating a high -resolution 
scatter x-ray image by extending said low-resolution scatter x-ray image D^ S | to all of said front low- 
energy detection locations through interpolation; and 

(k) calculating a high-resolution primary x-ray image D fHPh by subtracting said image D^,^ from said image 
Dfu, and calculating a high-resolution primary x-ray image by subtracting said image from said 
image D^; 

(I) whereby said images DfH Ph and Dfl_ Ph form a high-resolution, two-dimensional, dual-energy primary x-ray 
image pair D mPh , Dfl_p h of said subject after said scatter x-rays have been substantially eliminated, said image 
pair having a resolution substantially equal to the highest resolution availablef rom said front detector assembly. 
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cells not on projection lines, yielding 
scatter image D (HSh (x,y), and interpolate 
D fLSJ (x(i),y(j)) for front detector cells not 
on projection lines, yielding 
scatter image D 1LSh (x,y) 


Subtract D fHSh (x r y) from D fHh (x,y) to yield 
image D 1HPh (x,y) and subtract D fLSh {x,y) from 
D fLh( x > v ) t0 vie,d ima 9e D 1LPh (x f y) 


FIG. 12 
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